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Abstract 

Over the last decade several competing models of high-temperature super- 
conductivity were proposed, most of them with short-range interactions. We 
review a more realistic model with strong on-site repulsive correlations, the 
Coulomb and strong finite-range electron-phonon interactions. Bipolarons in 
the model exist in the itinerant Bloch states at temperatures below about 
half of the characteristic phonon frequency. Depending on the ratio of the 
inter-site Coulomb repulsion and the polaron level shift the ground state of 
the model is a polaronic Fermi (or Luttinger) liquid, bipolaronic high-Tc su- 
perconductor, or charge-segregated insulator for the strong, intermediate, and 
weak Coulomb repulsion, respectively. Two particular lattices are analysed 
in detail: a chain with the finite range electron-phonon interaction and a 
zig-zag ladder. Charge carriers in the ladder are superlight mobile intersite 
bipolarons. They propagate coherently without emission or absorption of 
phonons with about the same mass as single polarons. The model describes 
key features of the cuprates, in particular their Tc values, different isotope ef- 
fects, normal state pseudogaps, and spectral functions measured in tunnelling 
and photoemission. 
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I. INTRODUCTION 



For although high-temperature superconductivity has not been yet targeted as Hhe shame 
and despair of theoretical physics\ - a label attributed to superconductivity during the first 
half-century after its discovery - the parlous state of current theoretical constructions has led 
to a current consensus that there is no consensus on the theory of high-T^ superconductivity 
]![]. Our view is that the extension of the BCS theory towards the strong interaction between 
electrons and ion vibrations describes naturally the phenomenon, and the high temperature 
superconductivity exists in the crossover region of the electron-lattice interaction strength 
from the BCS-like to bipolaronic superconductivity p|. Quite remarkably Bednorz and 
Miiller noted in their original publication |^ and subsequently in their Nobel Price lecture 

that in their ground-breaking search for high-T^ superconductivity, they were stimulated 
and guided by the polaron model. Their expectation 0] was that if 'an electron and a 
surrounding lattice distortion with a high effective mass can travel through the lattice as a 
whole, and a strong electron-lattice coupling exists an insulator could be turned into a high 
temperature superconductor" . Indeed there is now overwhelming experimental and 
theoretical [pt p!0| , pT|p!2| , pr3tp!4| evidence for an exceptionally strong electron-phonon (e-ph) 
interaction in high temperature superconductors. Thus the theory of HTSC must include 
both e-ph and electron-electron Coulomb interactions as it was suggested some time ago . 
Also one has to take into account that all oxides are highly polarizable ionic lattices. A low 
density of mobile carriers is unable to screen effectively the direct Coulomb electron-ion and 
electron-electron interactions. The layered structure of the cuprates reduces screening even 
further. Since the mobile carriers are confined to the copper-oxygen planes their interaction 
with c-axis polarized optical phonons cannot be screened, and it is particularly strong. The 
parameter-free estimate of the polaron binding energy, Ep due to the long-range Frohlich 



e-ph interaction puts it at about 0.5 eV or larger in the cuprates [15 



However bipolaronic states are much heavier than band electrons since they are 'dressed' 
by the same lattice deformation, which bounds two polarons in a pair As a result. 
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the superconducting critical temperature, Tc, being proportional to the inverse mass of a 
bipolaron, might be significantly reduced rather than enhanced compared with the weak- 
coupling BCS Tc. Because of this prejudice some objections have been raised with respect 
to the bipolaron theory of high-temperature superconductivity. In the present article we 
review a few recent models of small bipolarons with the long-range electron-phonon (e-ph) 
interaction capturing the essential physics of superconducting cuprates [l^. These studies 



prove that small bipolarons are perfectly mobile Bloch states in the cuprates which can 
explain their high Tc. 

The electron-phonon coupling constant A is about the ratio of the electron-phonon in- 
teraction energy Ep to the half bandwidth D in a rigid lattice. We expect that when the 
coupling is strong, A > 1, all electrons in the Bloch band are "dressed" by phonons because 
their kinetic energy (< D) is small compared with the potential energy due to a local lattice 
deformation caused by an electron. If phonon frequencies are very low, the local lattice de- 



formation traps the electron. This self-trapping phenomenon was predicted by Landau |]T7 . 
It has been studied in greater detail by Pekar and Frohlich, and later on the most advanced 
path-integral theory of polarons was developed by Feynman and Devreese and his school in 
the effective mass approximation, which leads to the so-called large polaron (for more detail 



see Ref . [|TOl ) . The large polaron propagates through the lattice like a free electron but with 
the enhanced effective mass. In the strong-coupling regime, A > 1, the finite bandwidth be- 
comes important, so that the effective mass approximation cannot be applied. The electron 
is called a small polaron in this regime. The self-trapping is never "complete" , that is any 
polaron can tunnel through the lattice coherently. Only in the extreme adiabatic limit, when 
the phonon frequencies tend to zero, the self-trapping is complete, and the polaron motion 
is no longer translationally continuous. The main features of the small polaron were under- 



stood by Tjablikov Yamashita and Kurosava |T9|, Sewell [^], Holstein and his school 



21 1, Lang and Firsov |23], and others and described in several review papers and textbooks 



ZB| , |2^ J9| JT(]|J25|J2^J27| . The exponential reduction of the bandwidth at large values of A is 



one of those features. The small polaron bandwidth decreases with increasing temperature 



up to a crossover region from the coherent small polaron tunneling to a thermally activated 
hopping. The crossover from the polaron Bloch states to the incoherent hopping takes place 
at temperatures T uj^/ 2 oi higher, where is the characteristic phonon frequency. Here 
we show that small bipolarons are also in the itinerant Bloch states in a wide temperature 
region. Moreover, they propagate with about the same mass as a single polaron in particular 
lattices such as perovskites. 



II. SMALL POLARON BAND 



The canonical approach to a small polaron problem is based on the displacement (Lang- 



Firsov) transformation of the electron-phonon Hamiltonian in the site (m) representa- 



tion for electrons allowing for the summation of all diagrams including the vertex corrections, 

H = Y,T{in-n)5,,4c,+ (1) 
Y^uJc^hi [ui{q)dci + H.c] + ^ ^ V;(m - n)ni% + ^u^{dl^d^ + 1/2). 
Here T(m) is the bare hopping integral, 

is the matrix element of the electron-phonon interaction, i = (m, s), j = (n, s'), s =ti, 
hi = c-Ci, Ci,dq are the electron (hole) and phonon operators, respectively, and N is the 
number of sites {h = c = ks = !)■ 

Following Lang and Firsov one can apply the canonical transformation to diago- 



nalise the Hamiltonian. The diagonalisation is exact if T(m) = (or A = oo): 

H = e^He-^, (2) 

where 

S = -^hi[ui{q)dq^- H.c] (3) 
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The electron operator transforms as 



Ci = CiXi, (4) 



Xi = exp 



^Mj(q)(iq - H.c. 



and the phonon one as 



C^q = C^q - I]^i<(q)- (5) 



It follows from Eq.(5) that the Lang-Firsov canonical transformation shifts ions to new 
equilibrium positions. In a more general sense it changes the boson vacuum. As a result, 

ij i 

I]^q(4^q + 1/2) + \Y.'"i3^i^3^ (6) 



where 



CTij = T(m - n)6ss' exp ^Y^[uj{q) - Mi(q)](iq - H.c^ 
is the renormalised hopping integral depending on the phonon variables, and 



(7) 



Vij = Vcim - n) - ^ |7(q)|2wqCos[q • (m - n)] (8) 

q 

is the the interaction of polarons owing to the Coulomb repulsion ^0(111 — n) and to the local 
lattice deformation (the second term). 

In an extreme strong coupling limit A — > 00 one can neglect the hopping term of the 
transformed Hamiltonian. The rest has analytically determined eigenstates and eigenvalues. 
The eigenstates \N) — \ni, riq) are classified with the polaron 71^,3 and phonon occupation 
numbers and the energy levels are 

E = -Ep XI + I H ^ij^i^j + J2 ^q("'q + 1/2) (9) 

i ^ i^j q 

where = 0, 1 and riq — 0, 1, 2, 3, ....00. 
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Hence, the Hamiltonian Eq.(l) in zero order with respect to the hopping describes lo- 
cahsed polarons and independent phonons which are vibrations of ions relative to new equi- 
hbrium positions depending on the polaron occupation numbers. The phonon frequencies 
remain unchanged in this hmit. The middle of the electronic band falls by the polaronic 
level shift Ep as a result of a potential well created by the lattice deformation, 

^p=^El7(q)lV (10) 

With the finite hopping term polarons tunnel in a narrow band owing to the degeneracy 
of the zero order Hamiltonian with respect to the site position of a single polaron in a regular 
lattice. To see it one can apply the perturbation theory using 1/ A as a small parameter with 

A^^, (11) 

zT{a) 

where z is the coordination lattice number and T{a) is the nearest-neighbour hopping inte- 
gral, so that D Ri zT(a). The proper (Bloch) set of N degenerate zero order eigenstates of 
the lowest energy level {—Ep) of the unperturbed Hamiltonian is 

|k,0) = ^E4exp(^k-m)|0), (12) 

|0) is the vacuum. By applying the textbook perturbation theory one readily calculates 
the lowest energy levels of the polaron in a crystal. Up to the second order in the hopping 
integral the result is 

E(k) = -Ep + ek (13) 

V- l(k,0|E^,J■^^J^clcj|k^nq)P 

where |k',nq) are exited states of the unperturbed Hamiltonian with one electron and at 
least one real phonon. The second term in Eq.(13), which is linear with respect to the bare 
hopping, determines the small polaron band dispersion as 

ek = $]t(m)e-^'('")exp(-ik-m), (14) 



with the band-narrowing factor (at zero temperature) 



/(m) = ^ E l7(q)l'[l - cos(q • m)] (15) 



2N ^ 



The third term in Eq. (13), quadratic in T(a), yields a negative "k- independent correction 
to the polaron level shift of the order of 1/A^. 

III. SMALL BIPOLARON BAND 

The attractive energy of two small polarons, 2D{\ — yUc) is generally larger than the 
polaron bandwidth in the strong-coupling regime, 

A - Ate > ^' (16) 

because the bandwidth narrowing factor Z' is small 

^/ ^ Em T(m)e-g'('") exp(-zk ■ m) ^ ^ 
Em ^(m) exp(-zk ■ m) 

where /Zc is the Coulomb pseudopotential. When this condition is fulfilled, small bipolarons 
are not overlapped. Hence the polaronic Fermi liquid transforms into a Bose liquid of double- 
charged carriers. Here we encounter a novel electronic state of matter, a charged Bose liquid 
PJlI], qualitatively different from the normal Fermi- liquid and from the BCS superfluid. 



A. Onsite bipolaron band 

The small parameter, Z' j[\ — /ic) ^ 1, allows for a consistent treatment of bipolaronic 
systems p|j28[]. Under this condition the hopping term in the transformed Hamiltonian H is 
a small perturbation of the ground state of immobile bipolarons and free phonons, 

H = Ho + Hpert, (17) 

where 

Ho = IT. Vij4c]cjC, + uj^[d\d^ + 1/2] (18) 
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and 

Hpert = ^d-ijclcj (19) 



Let us first discuss dynamics of onsite bipolarons, wliicli are tlie ground state of the system 
with the Holstein-type non- dispersive e-ph interaction. The onsite bipolaron is formed if 

2Ep > U, (20) 

where U is the onsite Coulomb correlation energy (the so-called Hubbard U). The inter- 
site polaron-polaron interaction Eq.(8) is purely Coulomb repulsion because the phonon 
mediated attraction between two polarons on different sites is zero in the Holstein model. 
Two or more onsite bipolarons as well as three or more polarons cannot occupy the same 
site because of the Pauli exclusion principle. Hence, bipolarons repel single polarons and 
each other. Their binding energy, A = 2Ep — U, is normally larger than the polaron half- 
bandwidth, A ^ w = Z'D, so that there are no unbound polarons in the ground state. 
Hpert, Eq-(19), destroys bipolarons in the first order. Hence it has no diagonal matrix ele- 
ments. Then the bipolaron dynamics, including superconductivity, is described by the use 
of a new canonical transformation exp(S'2) ||2^, which eliminates the first order of Hpert, 

(C^^ ^ {f\^ij4cj\p) 

M- - Ep 

Here Ef^p and |/), \p) are the energy levels and the eigenstates of Hq. Neglecting the terms 
of higher orders than (w/A)'^ we obtain 

Wfr^ie'^He-'^^)^^,, (22) 
iHb)ff' ~ {Ho)ff' - E {f\^ii'ch'\p){p\^n'cyj'\f') X 

( \ \ \ 



\Ep-Ef Ep-EfJ 

S2 couples a localised onsite bipolaron and a state of two unbound polarons on different 
sites. The expression (22) determines the matrix elements of the transformed bipolaronic 
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Hamiltonian Hb in the subspace |/), |/') with no single (unbound) polarons. On the other 
hand the intermediate bra {p\ and ket \p) in Eq.(22) refer to configurations involving two 
unpaired polarons and any number of phonons. Hence we have 



where n{^^ are phonon occupation numbers (0, 1, 2, 3...oo). This equation is an explicit 
definition of the bipolaron binding energy A which takes into account the residual intersite 
repulsion between bipolarons and between two unpaired polarons. The lowest eigenstates of 
Hf, are in the subspace, which has only doubly occupied c|^|c|jjJO) or empty |0) sites. Onsite 
bipolaron tunnelling is a two-step transition. It takes place via a single polaron tunneling to 
a neighbouring site. The subsequent tunneUing of its "partner" to the same site restores the 
initial energy state of the system. There are no real phonons emitted or absorbed because 
the bipolaron band is narrow. Hence we can average Hh with respect to phonons. Replacing 
the energy denominators in the second term in Eq.(22) by the integrals with respect to time, 




(23) 




we obtain 



E E T(m - m')T(n - n') X 



(24) 




Here $iilin/(t) is a multiphonon correlator. 




(25) 



where 



Xi{t) = '[[exp[ui{q,t)dq - H.c], 



q 



and Mi(q,t) = ■Ui(q)e*'^''*. Xj(t) and Xj'(t) commute for any 7(q,i^) = 7(— q,i^) if m 7^ m'. 
Also Xj and Xj' commute, if n n', so that we can write 
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q 



Applying twice the identity 



yields 



Xj{t)Xiit)X]Xj, = He^'^'^e-^'^'^e 



f3'-dl-f3d^-\l3\y2 



X 



^[«,{q,t)-«,, (q,t)] [«* (q)-«*, (q)]/2-H.c. 



where 



(26) 
(27) 



(28) 



Finally using the average [22 



where = [exp(ti;q/T) — 1] ^ is the Bose- Einstein distribution function of phonons, we find 



^^P ^ l7(q)P^q(m, m',n,n') 



(29) 



cosh 




,2T 




sinh 


2T 





where 



Fq(m, m', n, n') = cos[q ■ (n' — m)] + cos[q ■ (n — m')] - 
cos[q ■ (n' — m')] — cos[q ■ (n — m)]. 



(30) 



Taking into account that there are only bipolarons in the subspace, where Hb operates, 
we finally rewrite the Hamiltonian in terms of the creation bj = cl^^cl^^ and annihilation 
bi = CmiCmi singlet pair operators as 



10 



A + - X! ^^^'Hm - m') 



1 



nm + 



(31) 



t(m - va!)h^j3^, + -v(m - m') 



There are no triplet pairs in the Holstein model, because the Pauli exclusion principle does 
not allow two electrons with the same spin occupy the same site. Here = bl^b^ is the 
bipolaron site-occupation operator, 



v{m — m') = 411(111 — m') + v^'^\m — m'^ 



(32) 



is the bipolaron-bipolaron interaction including the direct polaron-polaron interaction 
v{m — m') and a repulsive correction of the second order in T(m), 

,(2) 



v^^^ (m - m') = 2% / dte-'^*$™;^(t) 



(33) 

This additional repulsion appears because a virtual hop of one of two polarons of the pair 
is forbidden, if the neighbouring site is occupied by another pair. The bipolaron transfer 
integral is of the second order in T(m) 



t(m - m') = -2zT2(m - m') / dte-'^^^'^.it). 



(34) 

The bipolaronic Hamiltonian, Eq.(31) describes the low-energy physics of strongly coupled 
electrons and phonons. We use the explicit form of the multiphonon correlator, Eq.(29), to 
calculate t(m) and ^(^^(m). If the phonon frequency is dispersionless, cUq = ujq, we obtain 



'^'mm'(^) = e~^^ ^""^""'^ exp [-2g\m - m')e 
'^^^'{t) = e-^^'^'""'"') exp [2g\m - m')e-*^o* 

at T <^ cuq. Expanding the time dependent exponents in the Fourier series and calculating 
the integrals in Eqs.(34) and (33) yield p9| 



t(m) 



2T^(m) ^ [-2^^(m)]' 



(35) 



and 
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V 




2T2 



(m) 



g-2(;2{m) 



oo 



(36) 



A 



1=0 



/!(1 + U/A)' 



When A ^ c^o, we can keep the first term only with / = in the bipolaron hopping integral, 
Eq.(35). In this case the bipolaron half-bandwidth zt{a) is of the order of 2w^/(zA). How- 
ever, if the bipolaron binding energy is large, A ^ cuq, the bipolaron bandwidth dramatically 
decreases proportional to e~^^^ ^ 1 in the limit A — oo. However, this limit is not realistic 
because A = 2Ep — K < 2g'^UQ. In a more realistic regime, Uq < A < 2g'^ujQ, Eq.(35) yields 



On the contrary, the bipolaron-bipolaron repulsion, Eq.(36) has no small exponent in the 
limit A — s> oo, f*-^-* cx Z^^/A. Together with the direct Coulomb repulsion the second order 
f ensures stability of the bipolaronic liquid against clustering. 



Onsite bipolarons are very heavy for realistic values of the onsite attractive energy 2Ep 
and phonon frequencies. Indeed, to bind two polarons on a single site, 2Ep should overcome 
the onsite Coulomb energy, which is typically of the order of 1 eV or higher. Optical phonon 
frequencies are about 0.1-^0.2 eV in novel superconductors like oxides and doped fuUerenes. 
Therefore in the framework of the Holstein model, the mass enhancement exponent of onsite 
bipolarons in Eq.(37), is rather large (> exp(2£'p/co'o) > 150), so that onsite bipolarons could 
hardly account for high values of the superconducting critical temperature [Q. 

But the Holstein model is not a typical model. The Frohlich interaction with optical 
phonons, which is unscreened in polaronic systems, is much stronger. This longer-range 
interaction leads to a lighter polaron in the strong-coupling regime. Indeed, the polaron is 
heavy because it has to carry the lattice deformation with it, the same deformation that 
forms the polaron itself. Therefore, there exists a generic relation between the polaron 
stabilization energy. Bp, and the renormalization of its mass, m oc exp (7£'p/co'q), where 
the numeric coefficient 7 depends on the radius of the interaction. For a short-range e-ph 




(37) 



B. Intersite bipolaron band in a chain model 
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interaction, the entire lattice deformation disappears and then forms at another site, when 
the polaron moves between the nearest lattices sites. Therefore, 7 = 1 and polarons and on- 
site bipolarons are very heavy for the characteristic values of Ep and ujq. On the contrary, in 
case of a long-range interaction, only a fraction of the total deformation changes every time 
the polaron moves and 7 could be as small as 0.25 [^. Clearly, this results in a dramatic 
enlightening of the polaron since 7 enters the exponent. Thus the small polaron mass could 
be < 10 rrie where a Holstein-like estimate would yield a huge mass 10, 000 rrie. The lower 
mass has important consequences, because lighter polarons are more likely to remain mobile 
and less likely to trap on impurities. 

The bipolaron also becomes much lighter, if the e-ph interaction is longer-range. There 
are two reasons for lowering of its mass with increasing radius of the e-ph interaction. The 
first one is the same as in the case of a single polaron discussed above. The second reason 
is the possibility to form intersite bipolarons, which, in certain lattice structures, tunnell 



coherently already in the first-order in T(m) ||3^ , in contrast with onsite bipolarons, which 
tunnel only in the second-order, Eq.(35). 

To illustrate the essential dynamic properties of bipolarons formed by the longer-range 
e-ph interaction let us discuss a few simplified models. Following Bonca and Trugman 



we first consider a single bipolaron in the chain model of Ref . ^1| . One can further simplify 
the chain model by placing ions in the interstitial sites located between Wannier orbitals of 
one chain, and allowing for the e-ph interaction only with the nearest neighbours of another 
chain, as shown in Fig.l. The Coulomb interaction is represented by the on-site Hubbard U 
term. The model Hamiltonian is 

H = T{a) Y.i4+i,sCj-^ + ^-c-] + ^0 E 9{hj)njs{4 + d,) + (38) 

ujo J2[4^i + 1/2] + f^Xl ^jr^ji 



in the site representation for both electrons and phonons , where 
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and j are integers sorting the ions and the Wannier sites, respectively. This model is 
referred as the extended Holstein-Hubbard model (EHHM) [^]. We can view the EHHM 
as the simplest model with longer range than Holstein interaction. In comparison to the 
Frohlich interaction the EHHM lacks long-range tail in the e-ph interaction, but reveals the 
similar physical properties. In the momentum representation the model is a one-dimensional 
case of the generic Hamiltonian, Eq.(l), with 



7(g) = 9oV2{l + e^"'^) (39) 



and uj{q) = uq. Using Eqs.(lO), (14) and (8) we obtain 



^ TT i-7r/ 

for the polaron level shift. 



/7r/a 
dq[l + cosga] = 2(7qu;o (40) 
-7r/a 



g' = ^-^r dq[l - cos' qa] = g', (41) 

71 J-n/a 

for the mass enhancement exponent, and 



viO) = U - Agluo, (42) 



via 



/TT/a 
dq[l + cosga] cosga = —2gQUJo 
-n/a 



TT J-tt/o 

for the onsite and intersite polaron-polaron interactions, respectively. Hence the EHHM has 
the numerical coefficient 7 = 1/2, and the polaron mass 

m*EHP « exp (^^^ (43) 

scales as the square root of the small Holstein polaron mass, m*g^p a exp(i?p/c(Jo). In the 
case when U < 2gQUo, the onsite bipolaron has the lowest energy because |f (0)| > |f (a)|. In 
this regime the bipolaron binding energy is 

A = 4g',uJo - U. (44) 

Using expression (34) for the bipolaron hopping integral we obtain the bipolaron mass as 
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rriEHB oc exp ( ^ ) , (45) 



if A ^ cuq. It scales as {m*^fjp/mY , but occurs to be much smaller than the onsite bipolaron 
mass in the Holstein model, rnsuB °^ exp(4£^p/a;o), which scales as {m*gup/mY. In the 
opposite regime, when U > 2gQU!o, the intersite bipolaron has the lowest energy. Its binding 
energy 

A = 2gluo (46) 

does not depend on U. Different from the onsite singlet bipolaron, the intersite bipolaron 
has four spin states, one singlet S = and three triplet states, S = 1, with different z- 
components of the total spin, 5"^ = 0, ±1. In the chain model, Fig.l, the intersite bipolaron 
also tunnells only in the second order in T(a), when one of the electrons within the pair 
hops to the left (right) and then the other follows. This tunnelling involves the multiphonon 
correlation function Eq.(29), 

Hence the intersite bipolaron mass enhancement is 

m**„^ cx T-\a) exp (^) cx (47) 

in the infinite Hubbard U limit, C/ — > oo. We see that the intersite bipolaron in the chain 
model is lighter than the onsite bipolaron, but still remains much heavier than the polaron. 



C. Superlight intersite bipolarons in a ladder 

Any realistic theory of doped ionic insulators must include both the long-range Coulomb 
repulsion between carriers and the strong long-range electron-phonon interaction. Prom 
theoretical standpoint, the inclusion of the long-range Coulomb repulsion is critical in en- 
suring that the carriers would not form clusters. Indeed, in order to form stable bipolarons, 
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the el-ph interaction has to be strong enough to overcome the Coulomb repulsion at short 
distances. Since the el-ph interaction is long-range, there is a potential possibility for cluster- 
ing. The inclusion of the Coulomb repulsion Vc makes the clusters unstable. More precisely, 
there is a certain window of Vc/Ep inside which the clusters are unstable but bipolarons 
nonetheless form. In this parameter window bipolarons repel each other and propagate in a 
narrow band. At a weaker Coulomb interaction the system is a charge segregated insulator, 
and at a stronger Coulomb repulsion the system is the Fermi liquid, or the Luttinger liquid, 
if it is one-dimensional. 

Let us now apply a generic "Frohlich- Coulomb" Hamiltonian, which explicitly includes 
the infinite-range Coulomb and electron-phonon interactions, to a particular lattice structure 
33| . The implicitly present infinite Hubbard U prohibits double occupancy and removes 



the need to distinguish the fermionic spin. Introducing spinless fermion operators Cn and 
phonon operators dmu, the Hamiltonian is written as 

H=Y.T{n- n')cic^, + ^ V,{n - n')cicj^,c^, + (48) 

Here the e-ph and phonon terms are written in real space, which is more convenient in 
working with complex lattices, (7i/(m — n) is a dimensionless force acting between the electron 
on site n and the displacement of ion m, em-n = (m — n)/|m — n| is the unit vector in 
the direction from the electron m to the ion n, and is the polarization vector of the 
phonon branch u. Atomic orbitals of an ion adiabatically follow its motion. Therefore the 
electron does not interact with the displacement of the ion, whose orbital it occupies, that 
is gu{0) = 0. 

In general, the many-body model Eq.(48) is of considerable complexity. However, we 
are interested in the limit of the strong el-ph interaction. In this case, the kinetic energy is 
a perturbation and the model can be grossly simplified using the canonical transformation 
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of Section 2, which has the following form in the Wannier representation for electrons and 
phonons, 

m^n,i/ 

The transformed Hamiltonian is 

H = e^He-^ = Yl ^n„'cl.c„. + u;o E (^L^^m. + (49) 

n^n' ma ^ 

E f (n - n')cj,c„c|^,c„/ - -Ep ^ cj^Cn- 

The last term describes the energy which polarons gain due to el-ph interaction. Ep is the 
familiar polaron level shift 

Ep = u;Y9l{T^-i^){eu-e^^nY, (50) 

which is independent of n. The third term on the right-hand side in Eq.(49) is the polaron- 
polaron interaction: 

^;(n - n') = K(n - n') - 14,^(11 - n'), (51) 

where 

Vph(n - n') = 2ujQ ^ g^{m - n)5r^(m - n') x 

(g^ ■ Gm—n) i^^u ' ^m— n')' 

The phonon-induced interaction Vph is due to displacements of common ions by two electrons. 
Finally, the transformed hopping operator dnn' in the first term in Eq.(49) is given by 



<3'nn' = T{n — n') exp 



Y [g^T^ - n) (e,. • em-n) (52) 



- 5^^(111 - n')(e^ • e^-n')] - d^u) 

This term is a perturbation at large A. Here we consider a particular lattice structure 
(ladder), where bipolarons tunnell already in the first order in T(n), so that iTnn' can be 
averaged over phonons. When T <ujq the result is 
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t{n - n') = = T(n - n') exp[-/(n - n')], 



(53) 



g^{n-n') = ^ ^^(m - n)(e^ • 

[gu{ni - n){e„ ■ em-n) - gu{m - n'){e„ ■ Sm-n')] ■ 

By comparing Eqs.(53) and Eq.(51), the mass renormalization exponent can be expressed 
via E„ and V„h as follows 



2/ /\ 1 

g (n - n ) = — 



E,-lv,,{n-n') . (54) 



2 

Now phonons are "integrated out" and the polaronic Hamiltonian is 

Hp = Hq + Hpf,rt, (55) 

Ho = -EpY^cicn + J2 '^(^~ n')clcncl,Cn', 
n n^n' 

Hpert = J2 n')clcn'. 
n^n' 

When Vph exceeds Vc the full interaction becomes negative and polarons form pairs. 

The overall sign and magnitude of the interaction is given by the lattice sum in Eq.(51) 
evaluation of which is elementary. Notice also, that according to Eq.(54) an attractive 
interaction reduces the polaron mass (and consequently bipolaron mass), while repulsive 
interaction enhances the mass. Thus, the long-range character of the el-ph interaction 
serves the double purpose. Firstly, it generates additional inter-polaron attraction. This 
additional attraction helps overcome the direct Coulomb repulsion between the polarons. 
Secondly, the Prohlich interaction makes the bipolarons lighter. 

The many-particle ground state of Ho depends on the sign of the polaron-polaron inter- 
action, the carrier density, and the lattice geometry. Here we consider the zig-zag ladder. 
Fig. 2a, assuming that all sites are isotropic two-dimensional harmonic oscillators. For sim- 
plicity, we also adopt the nearest-neighbour approximation for both interactions, 5'i/(l) = g, 



Vc{n) = Vc, and for the hopping integrals, T{m) — Tjvat ior I — n — m — a, and zero 
otherwise. Hereafter we set the lattice period a — 1. There are four nearest neighbours in 
the ladder, z = 4. Then the one-particle polaronic Hamiltonian takes the form 

Hp^-EpJ2iciCn+pipn)+ (56) 
n 

Yl[t'(cl+lCn+pi+lPn) + t(picn + pl-iCn) + H.c], 
n 

where c„ and p„ are polaron annihilation operators on the lower and upper sites of the 
ladder, respectively. Fig. 2b. Using Eqs.(50) and (53) we find 

= Ag^coo, (57) 
7E„ 



The Fourier transform of Eq.(56) yields 

Hp = Y.{2t' cos k - Ep) {clck + plpk) + (58) 

k 

tY:[(l + e'')plc, + H.c.]. 

k 

A linear transformation of Ck and pk diagonalises the Hamiltonian, so that the one-particle 

energy spectrum Ei{k) is found from 

2t' cos k-Ep-Ei {k) t{l + 

t{l + e-'^) 2t'cosk- Ep - Ei{k) 

There are two overlapping polaronic bands. 



det 



0. (59) 



Ei{k) = -Ep + 2t'cos{k) ± 2tcos{k/2) 

with the effective mass m* = 2/|4t' ± t| near their edges. 

Let us now place two polarons on the ladder. The nearest neighbour interaction, Eq.(51) 
is V — Vc — Ep/2, if two polarons are on the different sides of the ladder, and v — Vc — 
Ep/i, if both polarons are on the same side. The attractive interaction is provided via 
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the displacement of the lattice sites, which are the common nearest neighbours to both 
polarons. There are two such nearest neighbours for the intersite bipolaron of the type 
A or B, Fig.2c, but there is only one common nearest neighbour for the bipolaron C, 
Fig. 2d. When Vc > Ep/2, there are no bound states and the multi-polaron system is a 
one-dimensional Luttinger liquid. However, when < Ep/2 and consequently v < 0, the 
two polarons are bound into an intersite bipolaron of the type A or B. 

It is quite remarkable that the bipolaron tunnelling in the ladder appears already in the 
first order with respect to a single-electron tunneUing. This case is different from both onsite 
bipolarons and from the intersite chain bipolarons discussed above, where the bipolaron tun- 
nelling was of the second order in T{a) . Indeed, the lowest energy degenerate configurations 
A and B are degenerate. They are coupled by i^pert- Neglecting all higher-energy configura- 
tions, we can project the Hamiltonian onto the subspace containing A, B, and empty sites. 
The result of such a projection is a bipolaronic Hamiltonian 

H, = (k - ^Ep) ^[AlA^ + BlB^] - t'^iBlA^ + B^i^n + ^-c], (60) 

where An — CnPn and Bn — PnCn+i are intersite bipolaron annihilation operators, and the 
bipolaron-bipolaron interaction is dropped (see below). Its Fourier transform yields two 
bipolaron bands, 

E2{k) = Vc--Ep±2t'cos{k/2). (61) 
with a combined width 4|t'|. The bipolaron binding energy in zero order with respect to t,t' 



IS 

A = 2E^{0) - E2i0) ^^-Vc. (62) 
The bipolaron mass near the bottom of the lowest band, m** = 2/t', is 

En 



m** =4:171* 1 + 0.25 exp ( ) . (63) 

The numerical coefficient 1/8 ensures that m** remains of the order of m* even at large 
Ep. This fact combines with a weaker renormalization of m*, Eq.(57), providing a superlight 
bipolaron. 

20 



In models with strong intersite attraction there is a possibihty of clasterization. Similar to 
the two-particle case above, the lowest energy of n polarons placed on the nearest neighbours 
of the ladder is found as 

6?2 — 1 

E„ = (2n-3)K (64) 

for any n > 3. There are no resonating states for n-polaron nearest neighbour configuration 
if n > 3. Therefore there is no first-order kinetic energy contribution to their energy. En 
should be compared with the energy Ei + {n — l)E2/2 of far separated {n — l)/2 bipolarons 
and a single polaron for odd n > 3, or with the energy of far separated n bipolarons for even 
n > 4. "Odd" clusters are stable if 

Tl 

and "even" clusters are stable if 

Tl — 1 

As a result we find that bipolarons repel each other and single polarons at > f-Ep. If Vc 
is less than ^Ep then immobile bound clusters of three and more polarons could form. One 
should notice that at distances much larger than the lattice constant the polaron-polaron 
interaction is always repulsive, and the formation of infinite clusters, stripes or strings is 
impossible Combining the condition of bipolaron formation and that of the instability 
of larger clusters we obtain a window of parameters 

^Ep < K < ^Ep, (67) 

where the ladder is a bipolaronic conductor. Outside the window the ladder is either charge 
segregated into finite-size clusters (small Vc) , or it is a liquid of repulsive polarons (large Vc) . 



IV. CONCLUSIONS 



It was analytically establshed a long time ago |TB|,pT],p2[ that small polarons are itinerant 
quasiparticles existing in the Bloch states at temperatures below the characteristic phonon 
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frequency for any strength of the electron-phonon couphng. Here we have reviewed more 
recent studies of small bipolarons in ionic lattices which show that the long-range Prohlich 
interaction leads to relatively light intersite small bipolarons with the atomic size of the 
wave function, large binding energy and a large size of the phonon cloud. They are Bloch 
waves with the effective mass, which is smaller by a few orders of magnitude than the mass 
of onsite bipolarons in the nondispersive Holstein model. As discussed in a few original 
papers, reviews and books the bipolaron theory describes Tc of many cuprates without any 
fitting parameters, their non Fermi-liquid normal state and the non-BCS superconducting 
state, and predicts single-particle spectral properties compatible with the tunnelling and 
ARPES spectroscopies of the cuprates. 

The author greatly appreciates stimulating discussions with A.R. Bishop, A.M. 
Bratkovsky, J.T. Devreese, D.M. Eagles, Yu.A. Firsov, L.P. Gor'kov, V.V. Kabanov, P.E. 
Kornilovitch, W.Y. Liang, K.A. Miiller, and S.A. Trugman. This work has been supported 
by the Leverhulme Trust (grant F/00261/H). 
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Figure captions 

Fig.l. Simplified chain model with two electrons on the chain interacting with nearest- 
neighbour ions of another chain. Second-order intcrsitc bipolaron tunneUing is sliown by arrows. 

Fig. 2. One-dimensional zig-zag ladder, (a) Initial ladder with the bare hopping amplitude 
T(a). (b) Two types of polarons with their respective deformations, (c) Two degenerate bipolaron 
configurations A and B. (d) A different bipolaron configuration C, which energy is higher than 
that of A and B. 
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